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The thermal Green function, T>^ v (x), for free, massless 
gauge bosons is computed exactly in a variety of gauges (Feyn- 
man, covariant, Coulomb, and Landshoff-Rebhan) . At large 
temporal separations it falls exponentially. At large spatial 
separations it falls like T/r. In contrast, the zero-temperature 
propagator falls quadratically in both regimes, being propor- 
tional to l/(t 2 - r 2 ). 
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I. INTRODUCTION 



Even though for fixed t this falls like 1/r 2 as r — > oc, it 
does contain the correct Coulomb potential, which comes 
from the light-cone singularity 5(t ± r)/r. At large time- 
like separations the 1/t 2 behavior causes the electron 
propagator to have a branch point at the electron mass 
shell instead of a pole. The 1/x 2 behavior often leads 
to processes that are divergent both in the ultraviolet 
(x — > 0) and in the infrared (x — > oc). 



A further consequence of Eq. (1.1) is that at zero tem- 
perature the quantum effects of massless particles are 
quite complicated because the effects are as important in 
space-like directions as in time-like directions. 



Although quantum field theories are formulated in 
four-dimensional coordinate space, calculations are al- 
most always performed in four-dimensional momentum 
space. One major exception to the preference for mo- 
mentum space is non-equilibrium field theory. Without 
equilibrium there is no invariance under time-translation 
and often no invariance under spatial translation. The 
dynamics of the density operator g controls the evolu- 
tion. For gauge bosons the time-ordered propagator is 

V^{x,y) = -iTr( Q T(A»(x)A»(y))). 

In non-equilibrium the correlations depend not just on 
the separation x — y but also on the effective age of the 
system x+y. For situations in which the dynamical evo- 
lution depends rapidly on x — y and slowly on x + y, it is 
standard to Fourier transform from the separation x — y 
to the conjugate four-momentum K. The transformed 
propagator is then a function of K and of x + y. It has 
long been known |j]-|4} how to use the Schwinger-Dyson 
equation for the non-equilibrium propagator to extract 
a kinetic equation for the non-equilibrium distribution 
function f(K) that depends on x+y . This approach has 
been extended to modern field theories containing gauge 
bosons particularly for QCD |fj-|{J. 

In ordinary vacuum field theory, or zero-temperature 
field theory, the density operator is g = |0)(0|. Transla- 
tion invariance is automatic and calculations are almost 
always performed in momentum space. In momentum 
space the time-ordered propagator for free gauge bosons 
in the Feynman gauge is (K) = — / \K 2 +irj) . This 
can easily be Fourier transformed to give the coordinate- 
space propagator 



^11 ( x )\t=0 



4tt 2 



le 



(1.1) 



Thermal Equilibrium 

For a field theory in an equilibrium heat bath that is 
uniform in space and constant in time, there is a constant 
temperature T throughout. The density operator is given 
by g = e~ H / T / r Tr(e~ H / T ). The time-ordered propagator 
has the structure ITol O 



V^{x) = 6{t)V^{x) + 



{-t)V^{x). 



(1.2) 



The thermal Wightman functions is defined by 

and P<(x) = £>>(— x). Appendix A summarizes how all 
other propagators (contour-ordered, retarded, advanced) 
can be expressed in terms of D > (x) . 

The thermal Wightman function satisfies two impor- 
tant conditions: Eq. (1.4) and (1.5) below. Both con- 
ditions come from the Heisenberg relation A^(t, r) = 
exp(iiit)A M (0, r) exp(— iHt). This implies that for com- 
plex time the field satisfies A' i (t,f) t = A*(t*,r). Con- 
sequently the Wightman function enjoys the reflection 
property 



[v^(x)]* = -v^(- 



(1.4) 



Furthermore, by inserting a complete set of energ y eig en- 
states between the two field operators in Eq. (|l.3|), it 
is simple to show that at complex time 2?>"(i, a?) is an- 
alytic in the strip —(3 < Imt < 0. A similar argu- 
ment shows that the Wightman function must satisfy the 
Kubo-Martin-Schwinger relation |13] 
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In momentum space the Wightman function is given 
by in terms of the momentum-space spectral function 
p^{K) by 



(1.6) 



which is derived in Appendix A. The free spectral func- 
tion is independent of temperature and only has support 
at k = ±|fc|. Thus only the on-shell particles contribute 
to the free Wightman function. The coordinate-space 
Wightman function is the Fourier transform 
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(1.7) 



As shown in Appendix A, the properties of the spectral 
function guarantee that X>> (x ) wi ll automatically satisfy 
the reflection condition in Eq. ( |l.4[ ) and the KMS relation 
in Eq. Q. 

The rema inde r of the paper will perform the computa- 
tion of Eq. (1.7) for free gauge bosons in various gauges. 
Sec II deals entirely with the Feynman gauge. The com- 
plete result as a f uncti on of arbitrary complex time is 
displayed in Eq. (2.14), but a simple r for m is the spe- 
cialization to real time shown in Eq. ( ^.15| ). Sec III per- 
forms the calculation in a general covariant gauge an d 
leads to the results summarized in Eqns. (3.6) and (3.7). 
Sec IV computes the Wightman function in the Coulomb 
gauge, which is a bit more subtle, and the results are 
given in Eqns. (O), (O), Q, and Oil). Sec V 



contains a comparison of propagators with Bose-Einstein, 
Fermi-Dirac, and Boltzman statistics. Appendix A pro- 
vides some general formulas. Appendix B performs the 
same calculation in the Landshoff-Rcbhan quantization 
scheme |l4|] in which the physical, transverse gauge fields 
are thermalized but not the longitudinal and time-like 
components. 



II. V^(x) IN FEYNMAN GAUGE 

It is easiest to compute the thermal Wightman func- 
tion in the Feynman gauge, in which the spectral function 
is proportional to the constant metric tensor g^ v : 

p^{K) = -g^27re(k )6(K 2 ). 

The thermal Wightman function has the form 

V^{x)=-grv > {x), (2.1) 

where 
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The scalar function T> > {x) is, in addition, the thermal 
Wightman function for a spinless field of zero mass. To 



compute the integral, first perform the integrals over ko 
and over the angles of k to obtain 



V>{x) 
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dkf k (r, t), 



(2.2) 



where the k dependence is contained in the function 



fk(t,r) 
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This function satisfies the conditions 

-fk(t,r)*=f k (-t*,r) 
f k (t-i0 >r )=f k (-t,r), 
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(2.3) 



which guarantee that T> > (x) wil l sat isfy th e re flection and 
KMS conditions given in Eqs. (O) and dO). 



A. Unregularized T> > (x) 

For any real value of t, except for t = ±r, the rapid 
oscillations of the integrand as k — > oo guarantee conver- 
gence. If t has a negative imaginary part then conver- 
gence is improved provided the negative imaginary part 
is smaller in magnitude than (3. Thus in the region 



(3 < Imi < 



t ^ ±r 



(2.4) 



the integral is convergent without any regularization. To 
perform the integration, it is convenient to use the fact 
that fk(t,r) is an even function of k in order to extend 
the integration range to the full k axis: 



2>>(a0 = 
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dkf k (r,t). 
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This can now be evaluated by contour integration over 
complex k. The integrand fk(t,r) is finite at k = 0. It 
has poles at k = ±i2nnT, for integer n > 1. 

How the contour should be closed depends on the rel- 
ative size of Ret and r. For definiteness, let Re£ > r. 
Then exp(±ifci) is the determining factor. For the term 
containing exp(— ikt) the contour should be closed in the 
lower-half of the complex k plane; for the term exp(+ikt), 
the contour should be closed above. Cauchy's theorem 
gives the result as a sum of the residues of the poles: 
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where the variables z± are defined by 



(2.6) 



This sum is only convergent in the region Re t > r. Per- 
forming the sum gives 
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2?>(i, r) 
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This result is valid for all complex t satisfying Eq. (2.4) 
and is analytic in this region. The poles at t = ±r and 
at t = ±r — i/3 (limit points of the open region) will be 
shifted slightly when the integration is regularized. It is 
convenient to express V > (x) in the alternate form 



P>(t,r) = 



-iT 
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coth(7rT(r+i)) + coth(7rT(r-f)) . (2 



Note that this is an even function of t. The regularized 
T> > (x) will include the domain t = ±r and will not be 
even under t — > — t. 

Deep time-like region: For |i| — r ^> 1/vrT, the asymp- 
totic behavior is 



D>(a: - — 
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Deep space-like region: For r — |t| 3> l/itT, the 
asymptotic behavior is not exponentially small. 
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(2.10) 



The dominant term T/r is present in all gauges and will 
be discussed further. 



B. Regularized T> > (x) 

The true Wightman function should be analytic in the 
closed region 

-f3<lmt<0, (2.11) 
which includes the points t = ±r and t = ±r — i(3 at 



which Eq. (2.7) has poles. The problem is that at these 



values of t the integrand fk{t,r) — > 1 without oscillation 
as k — > oo. The simplest way to r egula rize the integral 
in a way that will satisfy Eqs. (1.4) and (fLq) is to define 



-1 r°° 

T>>(x) = —j / dkf(t,r)e^\ (2.12) 
^ Jo 

where e is positive and real. Now the range of integration 
cannot be extended to negative k. However the integral 
can be evaluated using the relation flltf 
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= -4>(l + ai T) + i;(l + a 2 T), 



which holds whenever the real parts of ai and a 2 are 
positive. Here ip(z) = dh\Y(z)/dz . The result is 



2?>(x) = 



T 



ip[eT+iT(t-r)]-ip[eT+iT(t+r)] (2.13) 



-ip[l + eT-iT(t-r)]+ip[l + eT-iT(t+r)} 



This satisfies the reflection condition Eq. (1.4), which 
interchanges the first term with the second and inter- 
changes the thir d te rm with the fourth. It satisfies the 
KMS condition (1.5), which interchanges the first term 
with the third and the second term with the fourth. Us- 
ing ip( z ) = ^(1 + z ) ~ l/ z this can be written so as to 
isolate the zero-temperature contribution: 



T 
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ip[l+T(e-i(r+t))]-ip[l+T(e+i{r+t))] 



+ip[l +T(e-i(r-t))] - i/)[l+T(e+i(r-t))] 

(2.14) 
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This is the complete and gener al res ult. It is analytic in 
the closed region given in Eq. ( 2.11 ) . The nearest poles 
to this region are the zero-temperature poles just above 
the strip at t = ±r + ie and the the poles just below the 
strip at t = ±r — i((3 + e) from the first and fourth term. 

Real time: For real values of t one can safely put e = 
in all four tp functions since their poles are are at least 
a distance ±ij3 from the real t axis. When e = the psi 
functions have complex conjugate arguments and can be 
simplified using 



omp 




ip(l-iy) — ip(l + iy) — i7rcoth(7ry). 

The Wightman function is 
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coth[7rT(r + t) + coth[7rT(r - t) 

[8{r + t)-8[r-t)]. (2.15) 

This, of course, agrees with Eq. (2.8) when t ^ ±r. 
Note that under t — ► —t, the coth terms are symmetric 
but the Dirac delta terms are antisymmetric. The other 
Wightman function T> < (x) for real t is 



-iT 
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coth[?rT(r + t) + coth[7rT(r - t)} 
[5(r-t)-5(r + t)]. (2.16) 



Since Eqns. ( 2.15[ ) and ( p,16[ ) hold only for real t, one 
cannot pass from one to the other by the KMS relation. 



C. Time-ordered propagator 

The Wightman function satisfies the homogeneous 
equation □2?>(x) = 0. The time-ordered propagator is 

v n (x) = d(t)v > (x) + e(-t)v < ( x ), 

in which T> < (t,r) = T> > (—t ) r). This is a true Green 
function in that it satisfies the inhom ogeneous equation 
□ X>u(x) = ~5 4 (x). Using Eq. fl2.14Q gives 
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T> n (x)-- 



T 
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tp[l+T(e-i(r+t))]-ij[l+T(e+i(r+t))} 



+ 



+^[l+T(e-i(r-t))] - ij}[l+T{e+i{r-t))\ 
' (2-17) 
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This representation has the nice feature that the zero- 
temperature limit is isolated in the last term. The 
temperature-dependent terms are each annihilated by the 
d'Alembertian operator □ . 



D. Potential produced by a static charge 

In the limit r — > oo at fixed t, the Wightman func- 
tion and the time-ordered propagator have the behavior 
V{x) — > — iT/(Anr). This contribution suggests that at 
large distance the potential produced by a charge at rest 
would be temperature-dependent. This inference is in- 
correct as the following calculation demonstrates. 

Let J^{x) be a classical current. In the Feynman 
gauge, V^(x) = —g ,ll/ Vii(x) and so the classical vec- 
tor potential is 



A^(x) = -J d 4 x'D n ( 



x-x')J»{x'). 



For a point charge Q at rest J°(x') = QS 3 (f') and 
J(x') = 0. Thus the three-vector potential vanishes and 
the scalar potential requires integrating over the static 
charge density: 

/oo 
dt' V n (t-t',r) 
-oo 



It is convenient to use Eq. ( 2.17 ) for the Green function. 
At large complex t the combination of psi functions falls 
like 1/t 2 . Thus one can integrate over complex time by 
integrating over a contour C that is closed in the upper 
half-plane: 

A° cl {x) = -Q j> dt' V u (t - t', r) 



The first and fourth ip function in Eq. ( 2.17 ) are analytic 
in the upper-half of the complex t plane and make no 
contribution. The second and third ip functions have 
poles in the upper half-plane at t = ±r + i(N (3 + e) 

for N = 1, 2, 3 These poles all have the same residue 

(viz. —i/3) and their contributions to the potential cancel 
exactly. Thus the entire potential comes from the zero- 
temperature term: 



cA ' ~ 4ir 2 



dt' 



(t-t'f 



This is easily integrated and gives the usual Coulomb 
potential: 



(2.18) 



It is perhaps worth emphasizing that this not a large dis- 
tance approximation. A more difficult calculation, which 
does contain temperature-dependence, is the potential 
due to a point charge moving with constant velocity, 
where J°(x') = QS 3 {r' -vt'). 



III. 2?£ v (a:) IN CO VARIANT GAUGES 

The next case in which to compute the thermal Wight- 
man functions are the general covariant gauges. The 
tedious regularization performed in Sec II will not be 
attempted. The result given in Eq. (3.6) and (3.7) is 
therefore valid for real t not on the light cone. 

In a general covariant gauge the time-ordered propa- 
gator at zero temperature is 



The spectral function, extracted using Eq. (A7), is 



p^{K) = 27re(k )[-g^ + (1 - ^K^}S(K 2 ). 

When the thermal Wightman function T>^ (x) is ex- 
pressed in terms of the spectral function using Eq. (1.7) 
the result is 



V^(x) = -g^V > (x) + (l-0 



O 2 



dxudxv 



d>(x), (3.1) 



where T> > (x) with no superscripts is the Wightman func- 
tion from Sec II and the new function d> (x) is 



d > (x) 
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It will be important later that T>^(x) is not affected if a 
constant is added to the value of c?> (x) . The integration 
over ko and over the angles gives 



d>(x) 
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(3.2) 



where 
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The function gk{t,r) equals AT at k = 0. Consequently 
the integral in Eq. fl3.2j ) does not converge at k = 0. 
However the fact that dgk(t, r)/dt and dgk(t, r)/dr both 
vanish like k 2 as k — > guarantees that dd > (x) /dx u is 
finite, which is all that is necessary for Eq. (3.1). Since 
the behavior at k = is awkward, it is convenient to 



subtract a constant from Eq. (3.2) and redefine d > {x) as 
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d > (x) 



' dk u 
_[ 5fc (t, r )- 5fc (0,0)-^-^]. (3.3) 

The integrand now has no singularity at k = 0. A real 
parameter [i has been introduced so that the subtracted 
integration converges at k — oo for non-exceptional val- 
ues of t. Since it is not regulated as was done in Sec IIB, 
it does not converge for t = ±r o r fo r t = ±r — i/3. 

It is simple to integrate Eq. (|3.3|). The integrand is 
even in k and thus the range may be extended to — oo: 

,2 
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dk 
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[g k (t,r)-g k (0,0) 



k 2 + fi 2 - 



The integral can be evaluated by closing the contour in 
the complex k plane and using Cauchy's theorem. How 
the contour is closed depends upon the relative size of t 
and r. For example, if Ret > r the contour for the terms 
containing exp(— ikt) should be closed in the lower half- 
plane; those containing exp(ifct) sho uld be closed in the 
upper half-plane. The integrand of (3.3) is not singular 
at k = 0. The poles at k = ±i2irnT for n = 1,2,3.... 
give the result 
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where the constant terms independent of t and r are omit- 



ted and z± are the same as in Eq. (2.6). The summation 



displayed is only convergent if the magnitudes of z± are 
less than one: i.e. Ret > r. Performing the sum in this 
region gives 



d>(x) 



ln(l - z+) +hi(l - z-) 



(3.5) 
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This result can now be extended to Ret < r. Since it 
was not regularized, it is not valid precisely at z + — 1 or 
at z_ = 1. The result is not symmetric under t — > — i, 
but it can be rewritten as 



d>(x) 
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t + constant. 



(3.6) 



The term linear in i and the constant do not have second 
derivatives and thus do not affect the value of Df^(x). 
The function d>(x) is related to the scalar Wightman 
function by \3d>(x) = T> > (x). Therefore Eq. (3.1) can 
be summarized by 



-(x) = [-g^n+(l-0d^)d > (x). (3.7) 

The form of this in various limits will now be examined. 

Deep time-like region: Whe n |t| — r 3> 1/ttT the 
asymptotic behavior of Eq. (3J3) is 
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-2wT(\t\-r) , e -2wT(\t\+r) 



with an irrelevant term linear in t omitted. Consequently 
all components of (x) fall exponentially as was the 
case in Feynman gauge. 



Deep space-like region: If r — \t\ ^ 1/itT, then Eq. 
61) has the behavior 



iTr 

d >{x) -> -3— 
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omitting additive constants. The asymptotic behavior of 

,00 



V™(x) is 



V^{x)^§- r+ 0{e 



-2*Tr 



(3.8a) 



The leading term, being independent of ^, coincides with 
the Feynman gauge result. It was only this leading term 
that contributed to the calculation of the Coulomb po- 
tential in Sec II D. Of the remaining components, 2?> (x) 
is exponentially small: 



0(e 
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(3.8b) 



For the spatial components of the propagator, the term 
in d> (x) that is linear in r contributes a term from 



d 2 1 

r = -(5 13 -x l x 3 ). 



dx^dx j : 

The asymptotic form of the spatial propagator is 
IT 



V l i(x) -» — [-^+l)S^ + ^-l)x l x^\+0(e-^ lr ) 

(3.8c) 

Zero-temperature limit: At zero temperature, d>(x) 



is proportional to ln[r — t ] and Eq. (3.7 



gives 
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The combination x^x^Vft^ (x)\ T _ Q = £(£— 3)/87r 2 and 
thus vanishes in the Yennie gauge. 



IV. ©> J (x) IN COULOMB GAUGE 

In the Coulomb gauge the time-like component of the 
gauge potential is instantaneous and does not propagate 
in time. Consequently 

v^x) = S M 

47rr 

and there are no thermal corrections. For transverse 
gauge bosons the zero-temperature propagator is 



W \ 1 



k 2 J K 2 + ie' 
and the transverse spectral function is 

p lj (K) = 2TTe(k )S(K 2 )[S ij - Uti] 
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Using Eq. (1.7), the two-point function can be written 



D«(x) = S ij V > (x) + V 4 V J ff>(x), (4.1) 

where the first term, T> > (x), is the same as in Sec. II and 
the new function is 

H M--i ( * K e " K X < k ^ R2 ) 

By construction, T> > (x) = —V 2 H > (x). 

To compute if> (x) the first step is to integrate over kg 
and the angles to obtain 



, . -l r * , . , 
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where 
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For later reference note that H > (x) is analytic for com- 
plex time in the open region — (3 < Imt < 0. 

The function hk(t,r) equals 4iT at k = 0. Conse- 
quently the integral over k does not converge at k = 0. 
However dH > (x)/ dr i s convergent, which is all that is 
necessary for Eq. (4.1). To improve the k = behavior 
it is conve nien t to subtract an x-independent constant 
from Eq. (4.2) and redefine H > (x) as 



-1 



dk 



H>(x) = 8^ I w [Mi.r)-M0,0) 



k 2 



k 2 + /i 2 



where \i is an unimportant mass parameter. The integral 
is now convergent at k = and at k = oo. Since the 
integrand is even in fc, the range may be extended over 
negative k: 



H>(x) 
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(4.3) 



This will now be calculated explicitly for real t, without 
the regularization that was performed in Sec II. 

Time-like region: When t > r the contour for the 
e -ikt con tribution may may be closed in the lower half- 
plane and that for the e ikt contribution closed in the up- 
per half-plane. Although the full integrand has no poles 
at k = 0, each of these pieces has a double pole at k = 0. 
The more important contribution comes from the simple 
poles at k = ±i2irnT for n > 1. The result is 



H>{x) 



16ir 3 rT 



Li 2 (z_) - Li 2 (z+) 
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(4.4) 



The last two terms in come from the double pole at k = 
and do not contribute to dH > (x)/dr. The variables z± 
are given in Eq. ( [2.6| ) and Li 2 is the dilogarithm function, 
which for z < 1 has the series expansion |i~7| , |il| 



(4.5) 



When H > (x) is substituted into Eq. (4.1) the result is 

V${x) = {S 13 -x'x^Vyix) + (S l: > -3x l x 3 )E(x) (4.6) 
where the new function E{x) is 



E(x) 



^[ln(l-z_) + ln(l-z + )] 



167r 3 r 3 T 



;[Li 2 (z_)-Li 2 (2+)]. 



(4.7) 



In the deep time-like region, 7rT(|t| — r) 3> 1, both z + 
and Z- are exponentially small and so is E( x). Since 
T> > (x) vanishes exponentially as shown in Eq. ([^Jj) , this 
means that all terms in the Wightman function V l > (x) 
fall exponentially in the deep time-like region. 

Space-like region: When r > \t\ it is best to return 
to the defining integral Eq. (p~3|). In that integration, 
the contour for the e lkr contribution may may be closed 
in the upper half-plane and the contour for the e~ lkr 
contribution closed in the lower half-plane. As previously 
noted the full integrand has no poles at k = 0. However 
each of these pieces has a triple pole at k = 0. The simple 
poles at k = ±i2i:nT give infinite series: 



H>(x) 



167r 3 rT 
+ —\T(r 



Li 2 (— ) +Li 2 (z+) 



6rT J 



(4.8) 



The last term proportional to 1/rT is particularly un- 
usual. 

As mentioned previously, the original integral is ana- 
lytic in complex time provi ded — (3 < Imi < 0. However 
it is not obv ious that Eq. (4.8) is the analytic continua- 
tion of Eq. (4.4), and this provides an important check. 
In both results the arguments of the dilogarithm func- 



tions are smaller than one. For \z 
has the integral representation ll 



< 1 the dilogarithm 



tt 



Li 2 (z) = ds (In 
Jo 



s)- 



1 



(M<i) 



When the time is made complex by t — > t — ia then z± 



J2ttT<t 



z±\. Thus if < a < 0/2 then z± has a positive 
imaginary part. For \z\ > 1 the analytic continuation is 



U 2 {Z): 



f ds(lns) — —. r +in\n(z) ( 'f ' 

* - z I Al ' 



z\ > 1 
(z)>0 



Omission of the imaginary part, as is usually done in 
textbooks, would prevent the function from being ana- 
lytic. From the integral representation it follows that for 
\z\ > 1 and Arg(z) > the dilogarithm satisfies 



i„ , 2 n 



Li 2 (2) = y-~(mz)*-Li 2 (JJ +iirln{z). (4.9) 



G 



Using this relation one can analytically continue Eq. 
(4.4) from the time-like region where z_ < 1 to the space- 
like region where z_ > 1 and obtain Eq. ( Ell ) . Thus the 
two results agree despite their appearance. 

The Wightman funct ion in the spac e- lik e region results 
from substituting Eq. (4^8) into Eq. (4.1): 



where E(x) is now given by 



Zx i x j )E{x), (4.10) 



iT 

E ® = 8^ 



+ - 



87rr 3 

2 [Mi 



(Tt 2 



it ) 



l6ir 3 r 3 T 



[Li 2 (- 



- ln(l - 



(4.11) 



Naturally this is the analytic continuation of Eq. (4.7). 
It is easy to check that ViX>> (x) = 0. 

In the deep space-like region defined by wT(r— \t\) 3> 1, 
both l/z_ and z + are exponentially small. The asymp- 
totic behavior of the Wightman function is 



(x)^(^+i^)(^) 



(4.12) 



+ (S tj -3£ i x : >) 



i \ , 1 



(Tt 2 +it-— ) 



Despite its complicated appearance this still satisfies the 
transverse condition Vi2?> (x) = 0. The result is more 
complicated than in covariant gauges, which behaved as 
T/r with exponentially small corrections. Here the cor- 
rections are powers of r: the second line being order 1/r 3 
at fixed t but order 1/r if the ratio t/r is fixed as r— voo. 

Zero-temperature limit: It is straightforward to eval- 
uate the T = limit Coulomb-gauge Wightman function. 
As T — > the arguments of the dilogarithm functions in 
Eq. (4.7) approach unity: z — > 1. A useful expansion in 
this region is 0@ 

Li a (*) = ^-H0Mi-*)-£;^r-. (4.13) 

n— 1 

The Wightman function becomes 



T=0 



47r 2 (t 2 -r 2 ) 



-(5*^-34*^) 



8tt 2 





t + r- 








) 




U - r- 





(4.14) 



Naturally this agrees with the direct Fourier transform 
of the zero-temperature propagator. 



V. DISCUSSION 

The thermal Wightman function for free gauge bosons 
has been computed in various gauges. Knowing the 



Wightman function is the same as knowing the time- 
ordered propagators as shown in Appendix A. (The free 
retarded and advanced propagators are unchanged by 
the temperature since the free-field commutator is a c- 
number and thermal averaging does not change the c- 
number.) 

The rather surprising result is that the large-distance 
effects of massless gauge bosons are simpler at T > 
than in vacuum. In vacuum the free Green functions fall 
like 1 / (t 2 — r 2 ) . The thermal Green functions in covariant 
gauges and in the Coulomb gauge are exponentially small 
at large time-like separations. At large space-like sepa- 
rations the leading behavior is T/r, with exponentially 
small corrections in covariant gauges and power-law cor- 
rections in the Coulomb gauge. Appendix B computes 
the Wightman functions for the Landshoff-Rebhan ap- 
proach to thermalization |Q . 

Since the calculations presented all involve the Bose- 
Einstein distribution function, it is not apparent how 
much the asymptotic behavior depends on that particu- 
lar distribution function. The three examples below with 
different statistics will show that the results obtained in 
Sees II-IV are a consequence of quantum Bose statistics. 

Bose-Einstein Statistics: For Bose-Einstein statistics 
the simplest case is that of a massless scalar boson. The 
time-ordered propagator in momentum-space is 



Du(K) 



1 



2m8{K 2 ) 



K 2 



IT} 



1 



(5.1) 



The Fourier transform of this has been computed in the 
Feynman gauge discussion Sec II. For real times n ot o n 
the light cone, a useful representation is that of Eq. fl2.g| ): 



V>{x) 



-iT 



Airr 



coth(7rT(r+t)) + coth(7rT(r-t)) . (5.2) 



The asymptotic behavior as r — > oo can be understood 
by considering the propagator in the Matsubara formula- 
tion, which has discrete frequencies oj n = 2nirT because 
of the periodicity condition in imaginary time [|o| p^[ . 
As observed by Linde |h],|2(]] the n = modes produce 
the dominant behavior in the deep space-like region. The 
n = mode contribution is exactly —iT/^Anr). However 
the large time-like behavior does not come from a single 
mode in the Matsubara formalism, but requires summing 
a series and analytically continuing from complex ti me t o 
real time. The complete result is the same as Eq. (5.2). 

Fermi-Dirac Statistics: A convenient way to display 
the role of statistics is to change the distribution function 
in Eq. (5.1) from Bose-Einstein to Fermi-Dirac: 



Dn(K) = 



1 



2iri5(K 2 ) 



K 2 



IT] 



=/3|feol 



1 



(5.3) 



This differs from Eq. (5.1) in that the —1 in the denomi- 
nator has been changed to a +1 and also the overall sign 
of the second term has been changed. This is done so 
that the propagator for a massless spin 1/2 fcrmion is 
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the product of Eg (|5 3| ) with |HJ-|12|]. The Fourier 

transform of Eq. ( |5.3|) can be easily computed along the 
lines indicated in Sec II. For real times not on the light 
cone, the result is 



-iT 



Airr 



1 



1 



sinh(7rT(r + i)) sinh(7rT(r-i)) 



(5.4) 



At zero temperature this is the same as for Bose-Einstein 
statistics. At large time-like separations it falls expo- 
nentially as in the Bose-Einstein case. At large space- 
like s epa rations it also falls exponentially in contrast to 
Eq. (5.2). The space-like behavior can be understood in 
the Matsubara formalism, since the fermion frequencies 
Lu n = (2n + 1)ttT can never vanish. 

Classical Statistics: Classical Boltzman statistics 
gives a very different result. In the c lass ical limit the 
thermal distribution function in Eq. ( |5.l| ) is no longer 
singular at ko = 0: 



Dn(K) = 



1 



2md(K 2 ) 



K 2 



IT] 



=/3 1 fc 



(5.5) 



The Fourier transform of this gives 



T>>{x) 



—i 
4^2 



1 



1 



'■-(t-ie) 2 r 2 -{t-i(3) 2 r 2 -(t+i/3f 



At zero temperature this is the same as the Bose or Fermi 
case. However at T 7^ it is quite different. For large 
time-like separations it falls like 1/t 2 . For large space- 
like separations it also falls like 1/r 2 . Thus, for classical 
statistics the temperature does not substantially change 
the asymptotic behavior. 

Subsequent publications will explore the physical con- 
sequences of the space-time behavior. It should be pos- 
sible to understand the hard thermal loop approxima- 
tion |15[ ] directly in coordinate space. In that approx- 
imation the high temperature corrections come entirely 
from one-loop diagrams. The number of external lines 
determines the number of internal propagators. The im- 
portant effects come from thermal corrections on one 
internal propagator with all others kept at their zero- 
temperature value. A coordinate-space analysis should 
explain why only one loop diagrams are important and 
why only one propagtor in the loop enjoys thermal cor- 
rections. In addition, since non-equilibrium processes al- 
most demand a coordinant space treatment, it may be 
possible to deduce an extension of the hard loop approx- 
imation to non-equilibrium processes. 
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APPENDIX A: GENERAL STRUCTURE 

This appendix summarizes some standard properties 
of gauge boson propagators at finite temperature [^0|-fi2[ . 
The basic thermal Wightman function is 



2*T(*) = -*E e 



- E jn\A^x)A"(Q)\n) 
Tr[e-/3 ff ] 



(Al) 



It is customary to introduce a special notation T> < (x) for 
the function with inverted space-time arguments: 

V^{x)=V^{-x). (A2) 

In terms of the field operators this means that 

_ 0En (n\A»(0)A»(x)\n) 



V^(x) = -iJ2' 



(A3) 



Tr[e-0 H ] 

The retarded and advanced propagators are given by 

V%{x) = 0(t) \v^{x) - V^{x) 
V^(x) = 6(-t) \v^(x) - V^(x) 

For any parameter a in the range < a < f3, the four 
parts of the contour-ordered propagator are 

V^(x) = 9(t)V^(x) + 6(-t)V^{x) 

V^(x)=V^(t + ia,x) 

V^(x)=V^(t-ia,x) 

V^{x) = 6{t)V^(x) + e{-t)V^{x). 

With respect to real time, the first of these is time- 
ordered; the fourth is anti-time-ordered. 

The spectral function is the thermal average of the 
commutator 



-/3E„ 



(n\[A^{x),A-(0)]\n) 



Ti[e-P H ] 

Its Fourier transform p^ v {K) has the properties 

p» v (-K) = -p^(K) 
p^{K)* =p v »{K). 



The Fourier transform of Eq. ( |A1| ) is 

p^(K) 



1 



-f3k ■ 



(A4) 
(A5) 



(A6) 



This is the starting point in Eq. (1.6). The spectral 
function is most easily obtained from the propagator in 
momentum space by the relation 



r(K)=i[B%{K)-D>Z{lC)). 



(A7) 
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APPENDIX B: LANDSHOFF-REBHAN 
PROPAGATOR 

Landshoff and Rebhan have advocated heating only 
the two physical (i.e. spatially transverse) components of 
the gauge potential even in a general covariant gauge |Q . 
The procedure is to begin with the T — propagator in 
a general covariant gauge, 



-9 



K 2 



{K 2 



As explained in ||, the D m (K) and D°i(K) compo- 
nents are unaffected by temperature and therefore the 
corresponding Wightman functions in space-time may be 
read off from Eq. (3.9): 



> {X) = 8^ 



i (£ - 3)t 2 + (£ + l)r 



(;r k 



— ir]) 2 
tx J 
{x 2 — irj) 2 



(Bl) 
(B2) 



To separate the physical components of A, rearrange the 
spatial components of the T = propagator as 



D»(K) = 



5^ - k l W jk 2 (fcg - £k 2 )k l W jk 2 



K 2 



(K 2 + ie) 2 



The first term describes physical, transversely polarized, 
particles and is heated. The second term comes from 
gauge-dependent, longitudinally polarized particles and 
is not heated. From the previous results the finite- 
temperature two point function is 



V*{x) 



x) 



(B3) 



where the t emp er atur e-dependent Coulomb part is given 
by Eqs. (|4 .6 | ) , ( [t-7| ) in the time- like region and by 



Eqns. ( 4.10|) , ( 4.11 ) in the space-like region. The new 



temperature-independent contribution F> (x) is 
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- r 
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The Landshoff-Rebhan propagator clearly has a very 
complicated asympto tic behavior. In the de ep s pace-like 
region, r ->■ oo, Eqs. @ fal1 lik e V^ 2 , Eq. @ falls like 
1/r 3 , whereas the Coulomb part of Eq. ( fB3| ) still falls like 
1/r. In the deep time-like region, t — > oo, Eqns. (Bl), 



(B2), and (B3) are all dominated by the zero-temperature 
contributions and fall like 1 /t 2 rather than exponentially. 
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